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sociated with adiabatically generating a perturbed ground state are discussed as 
well. 
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1 Introduction 



Nonabelian gauge theories in general require the use of nonperturbative techniques as 
soon as one leaves the region of asymptotic freedom. Important examples of corre- 
sponding problems are the mechanism of confinement, probably related to the vacuum 
structure the investigation of hadron spectra and hadronic matter, the understand- 
ing of the deconfinement transition, the chiral symmetry restoration, the mechanism 
of hadronization in relativistic heavy-ion collisions, following the possible preformation 
of a quark-gluon plasma, and non-perturbative effects in a quark-gluon plasma above 
the critical temperature [Q. Besides numerically highly involved lattice calculations 
PI there are only preliminary attempts to achieve a nonperturbative description of 
nonabelian theories within the framework of Dyson-Schwinger equations [0-|lT[, dis- 



cretized light front quantization [0, variational calculations |T^-[|T^, and transport 
theories based on gauge-covariant Wigner functions with or without the restriction to 



an equal-time limit [1S]-|22 



In this article we follow a different approach along the line of n-point correlation 
dynamics which consists in a truncation of the hierarchy of equations of 

motion for equal-time Green functions in terms of connected Green functions up to 
order n. In relativistic field theory the method has previously been applied successfully 



to $^ theory in + 1, 1-1-1, and 2 + 1 spacetime dimensions [^, 29, |3^. The novel 
phenomenon in the application of ra-point correlation dynamics to Yang-Mills theory 
in temporal gauge is the requirement of Gauss law constraints on physical states. We 
apply 2-, 3-, and 4-point correlation dynamics to the simplest nontrivial Yang-Mills 
theory, i.e. SU(2) in 2 -|- 1 dimensions. The theory has the advantage of being not only 
superrenormalizable but finite [pl[] . 

This article is organized as follows: In Section ^ we discuss our procedure of quanti- 
zation in temporal gauge, which involves fixing the residual gauge freedom. In Section 
^ we present the derivation of the correlation dynamical equations of motion while in 
Section ^ we discuss the hierarchy of constraints for equal-time Green functions that 
arises from Gauss law. In Section |^ we present numerical results for the propagation 
of the system with given initial conditions for a variety of approximation schemes and 
investigate scaling, infrared and ultraviolet properties of the corresponding approxima- 
tions. In Section |^ we discuss the method of dynamically generating an approximate 
ground state configuration by the method of adiabatic switching. Section ^ contains a 
summary of our results. Appendix ^ contains a simple quantum mechanical example 
for our approach to fixing the residual gauge freedom, while in Appendix ^ we list the 
colour structure of the connected four-point functions within our approach. 
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2 Quantization in temporal gauge 



In this section we discuss the quantization of SU(N) Yang-Mills theories in temporal 
gauge {Aq = 0) and the corresponding free propagator. 

We first turn to the discussion of canonical quantization. The Lagrangian density 
of d + 1 dimensional SU(N) Yang-Mills theory is 

£ = -\f^,,F'^ , F;, = d,Al - d^Al + gr^^A^Al , (1) 

where a, 6, c are the colour indices in the adjoint representation of SU(N) and the greek 
indices run over the d + 1 spacetime components of the gauge field. Throughout the 
article, repeated indices of any kind will be summed over. The corresponding = 
gauge quantum Hamiltonian is 



H = jd'x {Ullimm + \F^,{S)Ft,{x)^ (2) 

and the canonical equal-time commutation relations are 

, t), n;(y, t)] = - y) , (3) 

where now the latin indices i, j only run over the d spatial components of the gauge 
field and 11^ denote the conjugate momenta. The Euler-Lagrange equations following 
from (|ip in the Aq = gauge for the spatial components of the gauge field A'} can 
be reproduced with (Q) and (0) as the Heisenberg equations for the operators and 
Uf. In contrast to that, the Euler-Lagrange equation for Aq, the Gauss law, is not 
contained in the dynamics generated by (0). It can, furthermore, not be realized as 
an operator identity, since this would contradict the commutation relations (0). The 
Gauss law thus has to be imposed as a constraint on physical states: 

G%x)\^,,) = , G%x) = (d^^Utix) + gr'^A^imm) ■ (4) 

One possibility to arrive at the above formulation of temporal gauge Yang-Mills theory 
is to directly deduce it from the Lagrangian (|l|), as we have just outlined, which involves 
right from the start to introduce a Hilbert space of wave functionals that depend on 
the spatial components of the gauge field, only. 

It is worthwhile to point out that, alternatively, one can start from a Lagrangian 
containing an additional (covariant gauge) gauge fixing term, and first introduce a 
Hilbert space of wave functionals depending on all -|- 1 spacetime components of the 
gauge field. The wave functionals then, in addition to Gauss law, have to fulfill the 
Lorentz condition, which can be used to obtain the above temporal gauge formulation 
by a separation of variables . 



The nontrivial aspects of the temporal gauge formulation now arise by the fact that 
the condition Aq = does not fix the gauge completely, and there is a residual gauge 
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symmetry consisting in the subgroup of all time independent gauge transformations 
(that tend to 1 at spatial infinity in the continuum or are periodic on a torus). The 
Gauss law operators C^^x) defined in (H) are the generators of this residual gauge 
symmetry, i.e. we have 

[G''{x),H] = (5) 

and for infinitesimal B"(x) (tending to zero at spatial infinity in the continuum or 
being periodic on a torus) 

c/-^[G]n,»(f)f/[e] = gr'^e^mK^) ■ (6) 

Condition (|D thus requires the physical states to be gauge invariant, i.e. 

U[Q]\^ph) = \'^ph) , (7) 

which also implies that the physical states are non-normalizable. 

The residual gauge freedom can now be dealt with in different ways. One possi- 
bility is to eliminate the remaining unphysical degrees of freedom and to formulate 
the theory entirely in terms of unconstrained variables. This can either be achieved 



by a change of variables |3^ or by applying unitary gauge fixing transformations 
p4| . A formulation in terms of unconstrained degrees of freedom as a starting point 
for nonperturbative calculations has the advantage, that the Gauss law (^) does not 
have to be observed as an additional constraint in an approximation scheme, since it 
is already implemented in the choice of variables. Furthermore, the wave functionals 
defined with respect to the unconstrained variables are then normalizable. However, 
such formulations have the disadvantage of leading to highly complicated (in general 
non-polynomial) expressions for the corresponding quantum Hamiltonians |Q. The 
application of correlation dynamics is, however, restricted to simple forms of the Hamil- 
tonian containing only low order polynomials in the fields. 

We will thus employ a different method for dealing with the residual gauge freedom 
©! ©) which consists in introducing a nontrivial integral measure for matrix elements 
between physical states obeying (^ by means of the usual Fadeev-Popov method (cf. 



P5 | and for a simple example Appendix We start by considering the matrix element 
of a gauge invariant operator Oinv (i.e., [Oi„^,G"(x)] = 0) between two physical states 
\%h) and \^ph): 

($ph|0,„,|vl>p,) = J Vf^[A] ^;,[A]Oi,^,,[A] , Ot^,,[A] ^ (A|0,„,|vl/,,) . (8) 

Since the integrand in (|^) is gauge invariant, in order to obtain a well defined result it 
is sufficient to choose the integral measure to contain a factor of the inverse of 



4 



the (infinite) volume of tlie residual gauge group, i.e. 



Vn[A] = i^J Vh[U]j VA , (9) 

so that the physical states can be chosen to be normalized in the following way: 

i^phl^ph) = J V^^[A] mi^[A]m,y,[A] = ( j Vh[U])-^ J VA %,[A]^,,[A] = 1 . (10) 

Here Vh[U] is the Haar measure for the integration over the residual gauge group and 
VA = l\g.i^^dA°;{x). The integral measure P/ifA], however, still leads to ill defined 
matrix elements for operators that are not gauge invariant. In order to retain the 
matrix elements of gauge invariant operators defined by and at the same time 

to have well defined matrix elements of gauge dependent operators, we now apply the 
Fadeev- Popov method to the integral (||). We first define a Fadeev- Popov determinant 
A[A] by 

1 = A[A] J Vh[U] 6[F[A^]] , (11) 

where F[A] =0 is some purely spatial gauge fixing condition and A^ denotes the field 
A transformed by the gauge transformation U. We then get in the usual manner 



{%H\0.n.\%h) = J mA] A[A] J Vh[U] mA'mH[APL%H[A] 
= J Vh[U] J Vf,[A^] A[A^]6[F[A^]]%,[A^]0fZ%,[A^] 

= jvA A[A]5[F[AWAA]0i^^AA] ■ (12) 

One could at this stage proceed to smear out the functional ^-distribution in the usual 
way, but this is not required for our further analysis. Eq. (0) defines an integral 
measure 

Vu[A] = VA A[A]6[F[A]] , (13) 

that yields the same results as Vfi[A] for the matrix elements of gauge invariant oper- 
ators and at the same time gives well defined finite results for the matrix elements of 
gauge dependent operators. 



The use of (|T3D for the definition of matrix elements in the space of physical states 
now has important consequences for the hermiticity properties of gauge invariant oper- 
ators such as the Hamiltonian (0). It is hermitean in the space of normalizable states. 



and retains this property also in the space of physical states [^, However, if it 
stands to the left of a gauge dependent operator (which leads out of the space of phys- 
ical states into the space of gauge dependent non-normalizable states), hermiticity will 
be lost (cf. Appendix^). 
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The standard example for the non-hermiticity of an operator, that is hermitean in 
the space of normahzable states, is the expectation value of the commutator A^{y)], 
where the (gauge-covariant) Gauss law operator loses its hermiticity ||35|, ^ |37| . 

In order to demonstrate the above mentioned non-hermiticity of the (gauge-invariant) 
Hamiltonian we consider the time evolution of the equal-time expectation value 
{^ph\dx^A°;{x)dyjA'j{y)\'^ph) in the abelian limit {g = 0). If we assume the system to 
be in an eigenstate of the Hamiltonian, 

H\^pj,) = E\^pj,) ^ H^^p^lA] = E^pj,[A] , (14) 
and work in the Schrodinger picture, we have 

'vv[A] (-iHHp^[A]y d,A^{x)dyA){y)^pn[A] 
+ / Vv[AWAA]d.A1{x)dy^A]{y) (-iHHp^[A]) = (15) 

(for simplicity, we use the same notation for classical fields and field operators). In 
contrast to that, the Heisenberg picture yields 

d,{%H\d,At{x)dy^A%y)\^p,,)\H.p. = {%H\][d.^At{x)dy^A'^iy), H]\%,) 

= {^p,\d.Mi^)dy,A%y)\<iIp,) + {^p,\d,Atix)dy^U%y)\^p,) 

= -id,^dy^6{x - y)6''%, = idf6{x - y)6^' , (16) 

where we have used (|^) and the abelian limit of (^). The discrepancy between ( P!3D 
and (|16D can be understood in terms of wave functionals, if one looks at the trans- 
formation from the Schrodinger picture to the Heisenberg picture at the time where 
the two pictures coincide: in order to get from one picture to the other, one has 
to partially integrate the functional derivatives appearing in the electric field energy 
I / d'^x [H"(x)n"(a;)]'^ = — ^ / d'^x sA^{x) ' 'which gives additional nontrivial terms 
when the derivatives act on the A-dependent terms appearing in the integral measure 



/ Vu[A] [/ d'zj^j^^%,[A]) d,At{x)dyA'^{y)%,[A] 



d^z j VA ^pnW X 
5 



^ 5Ai{z) 5Ai{z) [mmA]]d.AmdyA'AmvH[A] 



'd'^z J VA X 
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X 



^jvM\%.W ij d'^j^,jif^^^^mAmAA]) . (17) 



A different point of view leading to the same result is the following [^]: suppose instead 
of (1^) we would use (|]) and regularize the matrix elements of gauge dependent oper- 
ators by only integrating over a finite range of field amplitudes. The above mentioned 
additional nontrivial terms then arise as hermiticity-destroying surface contributions 
in the partial integration. 

In general, we note that the Heisenberg time evolution of gauge dependent quan- 
tities does not give the same result as the Schrodinger time evolution and yields time 
dependencies even if the system is in an equilibrium situation (e.g., an eigenstate of 
(H)). Also, in the Schrodinger picture the gauge fixing condition = is fulfilled 
at all times, while in the Heisenberg picture it is only fulfilled at the time of coinci- 
dence with the Schrodinger picture. However, if one considers the time evolution of 
gauge invariant quantities, both pictures give identical results (e.g., if one reconstructs 
an equilibrium expectation value of a gauge invariant quantity from gauge dependent 
quantities, all time dependencies in the Heisenberg picture must cancel out). Thus, 
since all physical obervables are gauge invariant, the Schrodinger and the Heisenberg 
picture are of course still equivalent. From a technical point of view, in order to generate 
equations of motion for Green functions from the Schrodinger picture, one would have 
to explicitly carry out the partial integration outlined above, which leads to nontrivial 
additional terms (containing the Fadeev-Popov (ghost) and gauge fixing contributions 
which arise due to the fixing of the residual spatial gauge freedom). 

One immediate consequence of the above considerations is the appearance of a term 
explicitly breaking time translation invariance in the free propagator as determined by 
the abelian limit of the Heisenberg equations of motion. If we choose the condition 
that fixes the residual gauge freedom to be the Coulomb gauge condition 

F[A] = d^Atix) = , (18) 

we obtain (in the continuum), if the Heisenberg and the Schrodinger picture coincide 
at time to; 



Xah I- Jdh. 



(27,) 



_}_^-i\k\\t-t'\ ( . _ 
]k\ V 



-o (l^-^'l + {t + t') -2t 







(19) 
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where T denotes the time ordering operator. A derivation of this result within the 
framework of canonical quantization can e.g. be found in 153, 53. We note in passing, 



that (|1^) at an intermediate time to 7^ =too might still not be a sufficient gauge fixing 
condition due to Gribov ambiguities. A detailed discussion of this issue, however, is 
beyond the scope of this article. 

The first authors to come up with a longitudinal propagator as contained in (p9|) 



were Carraciolo, Curci and Menotti |40], who obtained their result in a more heuristic 



way by first writing down the most general solution of the equation of motion. This 
solution, which is in addition required to be symmetric in t and t' and to fulfill the 
condition d'^{TAl i{x,t)Al/x\t')) = -dtdt>{TAl ^{x,t)A\ .{x',t')), reads 

{TAUS,t)AUS'.t')) = -'-{\t-t'\+a{t + t')+^)'- J ^,e^'''''^^. (20) 

In (|20| ) the subscript L stands for the longitudinal part of the gauge field. Now per- 
turbatively evaluating the Wilson loop up to order and comparing the final result 
to that obtained in Feynman and Coulomb gauge, the requirement of gauge invariance 
leads to the condition a = ±1, while 7 gives no contribution to the Wilson loop. The 
result a = — 1 is obtained in the above mentioned canonical derivation if instead of (^ 
one uses the condition 

i^phlG'^ix) = ^ G"^(f)^|vi>p,) = , (21) 



it 



which also leads to a consistent formulation of the theory. Note that G"(x) 7^ G"(x 
and that both (||) and (^) cannot be imposed simultaneously, since this would con- 
tradict d^). The authors in then conclude that contact with the Coulomb gauge 
can only be made at to = ice (or equivalently at one of the endpoints of a finite 
time strip). This is true if the Coulomb gauge condition requires every Green function 
containing a longitudinal gauge field Al^{x,to) to vanish. This, however, would imply 
that I7I has to be taken to infinity at some stage of the calculation. 

In order to shed some light on the last point we will conclude this section by briefiy 
discussing path integral quantization in temporal gauge. After carrying out the A^ = 
gauge fixing, one is left with a ghost-free path integral, that however still contains an 
integration over the residual gauge degrees of freedom. If one now fixes this residual 
gauge freedom via the Fadeev-Popov method with the Coulomb gauge condition at 
time to, i-e. 

d,A'}{x,to) = 0, (22) 



one arrives at 41 



{TAl,{x,t)Al/x',t')) = (It - t'\ - It - tol - It' - tol) ^ / ^d^''^'-''^ (23) 
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for the longitudinal part of the free propagator. Moreover, in the nonabelian theory 
the gauge fixing condition (|22|) leads to ghosts that live on a spacetime hypersurface 
of fixed time to. Eq. ( ^31) gives zero if to lies in between t and t' and it agrees with 
( PPD for a = —1 if t,t' > to and for a = +1 if t,t' < to. Agreement between 
and ([191) can thus be achieved only in the limit |to| oo. The disagreement between 
both results in case of a finite value of to, however, can easily be understood to arise 
from the difference in the way the Coulomb gauge condition is implemented in our 
canonical approach and in the path integral approach. The use of the integral measure 
(p!3[) induced by (0) for the definition of matrix elements between physical states only 
leads to vanishing Green functions if A1-{x, to) stands to the left of all other operators 
(for a simple example we refer the reader to Appendix while the implementation of 
(p!8D in the path integral formulation forces all Green functions containing A1 -(x, to) to 
vanish irrespective of the operator ordering (that arises from the time ordering being 
carried out). 

We therefore conclude that there is no a priori contradiction between (|T^) and (P3|), 
since both are in fact derived from inequivalent gauge fixing conditions. 

On one hand, we are now in a position to formulate simple Heisenberg dynamics 
employing the Hamiltonian (Q), on the other hand we have to deal with explicit time 
dependencies and the Gauss law (||) as an additional constraint. 



3 Equations of motion for connected equal-time Green 
functions 

We now turn to the derivation of a closed set of equations of motion for connected 
equal-time Green functions within the approximation scheme of correlation dynamics. 
The derivation of these equations up to the 4-point level has already been presented in 
|5[ for the full SU(N) QCD Lagrangian containing gauge bosons and massless fermions. 



however, in order to make this article self-contained we will repeat the main steps for 
the pure gauge boson sector in this section. We also point out that in the present 



article we no longer adopt the quantization scheme outlined in which is based 

on the "weak" version of Gauss law, i.e. requiring only ((^"(x)) = 0. 

The starting point for our derivation are the operator equations of motion that arise 
from the Heisenberg equation (if not explicitly stated otherwise, all time derivatives 
will from now on refer to the Heisenberg picture) 

idtO{t) = [0{t),H]. (24) 
With (|2|), d) and (|) one obtains 

dtA^ix) = U^ix) (25) 



9 



and 

d,Il^{x) = {v%,-d,^d,,)At{x) 

+2gr'^Alix)d,,Am + gr'^Amd.M^) - gr'^Alix)d.Alix) 
+g'r"'r'''Al{x)At{x)A<r{x) , (26) 

where all operators are taken at the same time t, which has been suppressed in our 
notation. In the same way one obtains the equations of motion for higher order equal- 
time operator products, e.g. 

d, {A^{x)A%y)) = {d,A^{x))A%y) + At{xMA%y)) 

= mx)A%y)+AKm%y), (27) 



dt {mx)A%y)) = {dmx))A';{y)+mxMA%y)) 
= - d,A,) Alix) 

+2gr''Al{x)d,^Af{x) + gr''A^{x)d.,Al{x) - g Al{x)d^Ai{S) 
+g'r'f''fAl{x)Al{x)A{{x)}A]{y) 

+n»(f)n;(y) , (28) 



= - d^^d,,) Atix) 

+2gr''Al{x)d,,Af{x) + gr''A^{x)d.,Al{x) - gr''Al{x)d,At{x) 
+g'r'f'^Alix)Alix)A{ix)}u%y) 

+Ut{x){{v%k-dy^dy,)Aliy) 
+2gf'^''Al{y)dy,Af{y)+gf^''Af{y)dy,Al{y)^ 

+g'f'''f'^Al{y)Al{y)A{{y)} , etc. (29) 

We thus arrive at an infinite coupled hierarchy of operator equations of motion, where 
the time derivative of each product of n field operators containing at least one conjugate 
field momentum couples to products of n, n + 1, and n + 2 field operators; the latter 
correspond to the free propagation and the gauge boson selfinteraction via the 3 point 
vertex and the 4 point vertex, respectively. 

Eqs. (p5D-(p9D and the corresponding equations for higher order operator products 
now straightforwardly induce an infinite hierarchy of equations of motion for full equal- 
time Green functions, the latter being defined as the expectation values of the above 
equal-time operator products. 
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The most general expression for the expectation value of an operator O is 



{0)=Tr{Op) , 

P = Y.Po'P^vh,^ ' = 1 ' Pa > Va (30) 

a a 

where P^Hp,^^ is the projection operator onto the physical state \^ph,a) and the set 
of states {\^ph,a)} is (without loss of generality) assumed to be an ortho-normalized 
basis system for the space of physical states, that diagonalizes the (non-equilibrium) 
statistical density operator p. The trace is only taken with respect to the space of 
physical states and each matrix element appearing in the trace is defined via the integral 
measure (|13|). The explicit form of (|30|) in terms of wave functionals then reads 

(O) = Ep- / Vi^[A]%^^MO^'^ph,M] ■ (31) 

Since we work in the Heisenberg picture, p is time independent and we have 

dtTr {0{t)p) = Tr {dtO{t)p) = Tr (^[0(t), H]p^ (32) 

for any operator 0{t). The equations of motion for full equal-time Green functions 
corresponding to (|^), (^1), etc., are therefore given by 

dt{Am) = m^)) , (33) 



+2gr'^{Al{x)d,^A'r{x)) + gr'^{A^{x)d.Xm - g {Al{x)d,^Al{x)) 
+g'r"'r''{Alix)Atix)Am) , etc. (34) 

The hierarchy consisting of (|33D, ( p^ ) and the corresponding higher order equations 
can be seen as the field theoretical analogue to the BBGKY (Born Bogoliubov Green 
Kirkwood Yvon) density matrix hierarchy known from nonrelativistic many body the- 
ory. The only new feature is the difference between the hierarchies arising from the 
Heisenberg and the Schrodinger picture, as it has already been discussed in Section ^ 

We note in passing that, in the Schrodinger picture, the von Neumann equation 
remains vahd also for the description of the time evolution of gauge dependent oper- 
ators, since for its derivation only the hermiticity of H in the space of physical states 
is needed (all manipulations in the derivation are carried out to the right of O). One 
could then be tempted to recover the Heisenberg picture result by using invariance 
of the trace under cyclic permutations. This, however, is incorrect, since the trace 
cannot be extended to run over a large enough space of states in order to allow for an 
application of the completeness relation as required for the proof of cyclic permutation 
invariance. 
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For practical calculations, the infinite coupled hierarchy of Heisenberg equations of 
motion for full equal-time Green functions has to be truncated. In our approach, this 



is done using the cluster expansion for ra-point Green functions [^, i.e. their decom- 
position into sums of products of connected Green functions. The cluster expansion is 
expected to exhibit good convergence properties with respect to the order of connected 
Green functions taken into account, if the system is in a single phase (localized wave 
functional) configuration. On the other hand, if the system is in a configuration with 
coexistence of different phases (e.g. in a groundstate wave functional that is highly 
delocalized due to tunneling), the higher order connected n-point functions are ex- 
pected to contribute to the full n-point functions in the same order of magnitude as 
the corresponding disconnected parts 



30 



The explicit form of the cluster expansion can be derived from the generating func- 
tionals of full and connected Green functions, Z[J, a] and W[J, a], given by 



Z[J,a] = Tr 
and Z[J, a] - 



P 



(35) 



where x = {x, t). We start with the cluster expansions of time ordered Green functions 
with different time arguments: 



(Atix)) = (A»(x)), , mx)) = (n»(x)). 



(36) 



{TArAxi)A^^ix2))= lim 



6 



.W[J,a] 



+ 



hm e^[''''^l 
6 



-W[J,a 



J,<^^ot5J^^^{x^)t6J^^'ix2) 
6 



jsj^,'{x2: 



W[J,a] 



\t6J^;{x,: 

(r<(xi)A-(x2))e+ (A-(X1))(A-(X2)) 



(37) 



where {■)c denotes the connected part of the expectation value. Analogously we obtain 

{TnrAxi)nt^ix2)) = {TnrAxiM^{x2)). + {nrAxi)){ut^ix2)) , 



(TA-(xi)A-(x2)A^3(x3)) = {TArAxMt!ix2)At^{xs))c 

+ {TAtl{xMt^{x2))Mt!{x,)) + {TAl^{xr)At!{x,))Mt^{x2)) 



+ (TA-(x2X(x3)).«(xi)) + {ArAxi)){At^ix2)){At:{xs)) , etc. 



(38) 
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The expressions for equal-time Green functions are obtained by taking the well defined 
equal-time limit which yields the desired operator ordering in the cluster expansions. 
We arrive at 

= {Anx)). , mx)) = mx)). , 

{At:iS,)A1^ix,)) = {ArAxMt^ix2))c+{At^ix^)){A1^ix,)) , 

{^t:{xi)At^{x2)) = {l^rA^Mt:i^2))c + {nt:ix^)){At^{x,)) , 
{iitAx,)iitAx,)) = (n^;(fi)nr,^(f2))c + (n-(fi))(n-(f2)) , 

{At:{xMtAx,)At:ixs)) = {ArAxMtAx2)AtAx,% 

+ (yl-(fi)^-,^(5^2))c(A-^(x3)) + {Al^{x^)A1^{x,))MtAS2)) 

+ {AtAx,)AtAx,MAt:{xi)) + {At:{x,)){AtAx2)){AtAxs)) , etc. , (39) 

where all time arguments are now equal and have been suppressed in our notation. In 
view of their length the cluster expansions of other Green functions required for our 
calculations are not explicitly given here, but e.g. can be found in |^ . 



We mention that instead of (pq), one could equivalently start with a generating 
functional that only depends on the source Jf (x) and then obtain Green functions 
containing conjugate momenta n"(x) by differentiating with respect to the time coor- 
dinate. An alternative for the derivation of (^) without the detour via time ordered 
Green functions is the use of an equal-time generating functional, where due to the 
absence of the time ordering operator one has to take care of the nonvanishing com- 
mutator between Af{x) and 11° (x) by rewriting the exponential containing the sources 
J and a by means of the Campbell-Baker-Hausdorff-formula [^ . 

Inserting the cluster expansions into the hierarchy equations of type (|33|), (^4|), etc., 
leads to an infinite coupled hierarchy of equations of motion for connected equal-time 
Green functions. Our approximation scheme now consists in truncating this hierarchy 
by neglecting all connected n-point functions with n > N, which leads to a closed set of 
coupled, nonlinear equations of motion for the connected equal-time Green functions up 
to the A^-point level. These equations are referred to as A^-point correlation dynamics. 
Since the time derivatives of the full A^-point functions couple to the full A^ + 2- 
point functions, all cluster expansions up to the A^ + 2-point level are required for the 
derivation of the correlation dynamical equations. 

In the truncation scheme with A^ = 1 we recover the classical Hamilton equations 
for temporal gauge Yang-Mills theory (the 1-point functions are simply the classical 
fields). In the N = 2 truncation scheme we recover the time dependent Hartree-Fock- 
Bogoliubov equations (this truncation scheme corresponds to a Gaussian trial state 
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and to a Gaussian form of the integral measure; in the present case this, however, 
contradicts gauge invariance). The = 3 truncation scheme is the lowest order ap- 
proximation to take into account the 3-point vertex if the system is in a (global) colour 
singlet configuration. 

The highest order approximation investigated in this article is = 4. For the 
straightforward but very lengthy derivation of the equations of 2, 3, and 4-point cor- 
relation dynamics we have developed computer algebra codes using the Mathematica 
programming language, which are also used to generate (highly optimized) Fortran rou- 
tines for the time integration of the equations within a finite set of plane wave states 
on a torus. The resulting equations of motion contain several hundred terms and are 
therefore not explicitly displayed in this article. The reader is, however, referred to our 
previous article which contains the explicit 4-point equations in coordinate space 
in a notation that is compactified using permutation operators. 



4 Hierarchy of Gauss law constraints 

In addition to the equations of motion following from the Heisenberg equation - dis- 
cussed in the preceding section - the equal-time Green functions have to fulfill the 
conditions that can be derived from Gauss law (H), which, as stated in Section ^, is 
not automatically contained in the Hamiltonian dynamics of the theory. The present 
section is devoted to the discussion of these Gauss law conditions. 
We first note that from (^) we immediately get 

OG'«(f)|^p,) = (40) 

for any product of field operators O and for any physical state \^ph)- For an arbitrary 
expectation value (^) we thus have 

{OG\x)) = , (41) 

which in principle already comprises the complete set of equations for equal-time Green 
functions that can be derived from Gauss law. It is important to point out that in ( PD 
G°'{x) stands to the right of all other operators and that the operator ordering can- 
not be reversed by invoking hermiticity arguments (the Gauss law operator, although 
hermitean in the space of normalizable states, loses this property in the present case 

If we replace O in (|4l|) by a product of n field operators, then (due to the structure 
of G"'{x)) we get a relation between a full n + 1- and a full n + 2-point equal-time 
function. Thus, as in the case of the Heisenberg equations of motion, we are left with 
an infinite coupled hierarchy of equations. In the simplest case, O = 1, we obtain a 
relation between 1- and 2-point functions: 

{G\x)) = d^Mm + gr'^Aiimm) = o . (42) 
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On the next level, with O = A'j{y) and O = Ilj{y), we get two different relations 
between 2- and 3-point functions: 



d^M^mtm + gniA^mKmm) = o , (43) 

d,AU%mtiS)) + 9r''%Il%y)A^{x)mx)) = . (44) 

Again one level higher we get three different relations between 3- and 4-point functions 
by choosing O = A'^{y)Al{z), O = U'^{y)Al{z), and O = U'^{y)Ul{z). We note that 
O = A'^{y)Ill{z) is redundant, since it can be obtained by commutation relations. 

d^M'mimti^)) + gr'^iA'^mimiimm = o , (45) 

d,^{U%y)Al{z)Ut{x)) + gr^W^mimlimm) = , (46) 

d.A^%mimtm + gr'^{n%mim-{mm) = . (47) 

Further relations are easily obtained for the higher orders. Within our temporal gauge 
formulation p^), (^3]), and the corresponding higher order equations are the 
Slavnov- Taylor identities for equal-time Green functions and have to be fulfilled in 
order to guarantee gauge invariance. 

If we consider the (Heisenberg picture) time derivative of one of the above Slavnov- 
Taylor identities, we have 

dtiOC^ix)) = {-[0,H]G''ix)) + {0-[G''ix),H]) = . (48) 



If O is an n-point operator, the first term in (48) vanishes due to the Slavnov- Taylor 
identities that relate the equal-time Green functions of order n + 1 and n + 2, those of 
order n + 2 and n -|- 3, and those of order n + Z and n + 4, since the time derivative 
of an n-point operator in general couples to n-, n -|- 1-, and n + 2-point operators. 
The second term vanishes due to (^). Thus, in principle, the hierarchy of Slavnov- 
Taylor identities can be fulfilled by an appropriate choice of initial conditions for the 
equal-time Heisenberg hierarchy. 

This last property, however, only holds if we consider an exact solution of the 
hierarchy of equations of motion; it breaks down for solutions within our approximation 
scheme. This comes about for two different reasons: 

a) Since for our numerical calculations we consider the theory on a torus and only 
take into account a finite number of plane wave states, identity (^ is violated, i.e. if H 
and G°'{x) are replaced by the corresponding finite momentum space representation, 
their commutator does not vanish any more. This is nothing but the well known 
fact that any simple momentum cutoff violates gauge invariance. In Section S we will 
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discuss how this can be partly remedied by considering different momentum cutoffs for 
the 2- and the 3-point functions. 

b) The truncation in terms of connected Green functions is not fully consistent with 
the constraint hierarchy following from Gauss law, so that even for a regularization 
scheme that preserves gauge invariance, not all Slavnov- Taylor identities would be 
conserved in time. 

In order to investigate the last point in more detail in the SU(2) case, we first note 
that the 2-, 3-, and 4-point correlation dynamical equations of motion - regardless of 
the finite momentum basis - conserve the following global colour singlet structure of 
connected equal-time Green functions (O denotes fields A and/or conjugate momenta 



n): 

mx)) = , (49) 

{Ot^{x,)Ol-{x,)), = S'^^'^mO,,{x,)OUx2)))c , (50) 

{0t^{x,)0l'{x2)0t!{xs))c = e'^^^-^^{0,,{x^)0,,{x2)0,,{xs)))c , (51) 



where we have /"''^ = e"'^'^ for colour SU(2). Here the double brackets ((■)) denote a 
colour independent function, depending only on the spatial coordinates and the vector 
indices of the corresponding field operators. The notation using operators without 
colour indices inside the double brackets is to be understood as purely symbolic and is 
only supposed to indicate the kind of field (A or 11) to which the spatial coordinates 
and vector indices belong (our notation should especially not be confused with that of 
matrix valued fields in the fundamental representation of SU(2)). While for the 2-point 
functions we have 

{{Ai,{xi)A,{x2)))c = {{Ai,{x2)Ai,{xi)))c , 

{{Ui,{x,)Ui,{x2)))c = ((n,,(f2)n,,(Xi)))e , (52) 

and 

{{Ai^{xi)Ui.^{x2)))c = {{Ili^{x2)Ai^{xi)))c + i5i^i25{xi - X2) , (53) 

the colour independent 3-point function, due to the antisymmetry of the SU(2) struc- 
ture constants, is antisymmetric with respect to the exchange of any pair of fields or 
momenta, i.e. we have 

((0,,(fi)0,,(f 2)0,3 (fa))) 

= -{{0,,ix2)0,,{x^)0,,{xs))) = -{{Oi,{x^)0,,{x3)0,,{x2))) . (54) 

The color structure of the connected 4-point functions is more involved and given 
in Appendix 0. Since our condition for fixing the residual gauge freedom (0) does not 
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prefer any direction in colour space, we always have such a global singlet configuration 
and can therefore make use of (pDD, (pT|), and ( |103| )-( ]TUB| ) for an explicit reduction 
of the number of degrees of freedom appearing in the correlation dynamical equations. 
This, as well as the proof of the conservation of the above colour structure by n-point 
correlation dynamics {n < A), is straightforward but lengthy and has been carried out 
by means of computer algebra. 

One immediate consequence of (|49| ) is that the full 2- and 3-point functions do not 
have disconnected parts and therefore equal the connected 2- and 3-point functions. 
Furthermore, (^21) is trivially conserved within correlation dynamics, since for all times 
we have 

(G«(f)) = d^Mi^)) + ge^'^'iA^imK^)) = ^?e'^'^5^^p,(f)n,(f))), = O . (55) 



Also, due to the global colour singlet structure, identities p3| ) and (0), relating 
full 2- and 3-point functions, directly become identities for connected Green functions. 
The Gauss law constraints of higher order, relating full n- and n + 1-point functions, 
in general translate into identities relating connected n- and n + 1-point functions and 
products of all lower order connected equal-time functions down to the connected 2- 
point function. Thus, within an n-point truncation scheme, all identities up to the 
one relating the full n — 1- and n-point functions remain unchanged, and the identi- 
ties relating the full n- and n + 1-point functions have to be truncated by neglecting 
the connected n + 1-point function; e.g. within a 3-point truncation scheme one has 
relations between connected 2- and 3-point functions and identities relating connected 
3-point functions to products of 2 connected 2-point functions. 

Since the conservation of a Gauss law constraint identity, as stated above, requires 
the conservation of the identities of the 2 next higher orders, a truncation in terms of 
connected Green functions up to a certain order in general leads to a non-conservation 
of (15), (|4^) , and the higher order constraints. 



A subset of conserved identities can nevertheless be obtained if, instead of multiply- 
ing G°'{x) from the left with an arbitrary operator product, one considers expectation 
values of products of several Gauss law operators. In general, a product of n Gauss law 
operators contains products of up to 2n fields or momenta, and its expectation value 
is thus conserved within a 2n-point truncation scheme as long as [G"'{x),H] = P7|. 



The conservation of these identities is of course destroyed as soon as we work within a 
finite momentum basis. 

It is worthwhile to mention, that within the framework of correlation dynamics for 
nonrelativistic fermion systems there exists a constraint, that is violated in a similar 
way by the n-point truncation scheme as the Gauss law in the present case of correlation 
dynamics for a nonabelian gauge theory. This constraint is the requirement that the 
system has good particle number, i.e. if Npermi = J2a (^a(^a is the number operator 
and Hpermi is the Hamiltonian, with [Hpermi, Npermi] = 0, and I'^Permi) is the state 
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describing the system in Fock space, we impose the constraint 



(^-^ Fermi {-^ Fermi )) l"^ Fermi) = , (56) 

from which we directly obtain 

(O {NFerm^ ' {NFerm^))) = (57) 

for any operator O. By inserting any possible product of creation and annihilation op- 
erators for O we obtain a constraint hierarchy similar to the hierarchy of Slavnov- Taylor 
identities as obtained above from Gauss law. This constraint hierarchy is nothing but 
the set of trace relations relating the n-body (2n-point) and the n + 1-body (2n + 2- 
point) density matrices of a pure iV-particle wavefunction, which is well known to be 
violated within the correlation dynamics approach. Just as in the case of Gauss law 
a subset of conserved identities within a 2n-point truncation can be obtained by con- 
sidering the expectation values of -/V^ermi; k < n (e.g. (Npermi) conserved within 
2-body (4-point) correlation dynamics). 

However, there is still a major qualitative difference between the Gauss law in Yang- 
Mills theory and the trace relations in fermionic many-body theory. While the Gauss 
law (^) is induced by a local gauge symmetry, condition (p6D is induced by a global 
symmetry. Therefore particle number conservation is not violated on the operator level, 
if the theory is considered within a finite single particle basis. 

While the violation of the trace relations was found to be only weak and without 
significant impact on the time evolution of the system at least for the application of 
correlation dynamics to light nuclei it is not a priori clear if this will also be the 



case with the violation of Gauss law in temporal gauge Yang-Mills theory and has to 
be studied explicitly. 



5 Numerical investigations for SU(2) in 2+1 di- 
mensions 

For our numerical calculations we consider SU(2) and expand the fields and conjugate 
momenta into plane wave states on a torus in 2 + 1 spacetime dimensions according to 

AKx) = E y^'Atik) , n^(f) = E y^^^Utik) , (58) 

k k 

where L is the size of the 2-dimensional box (of volume L^) and the momentum vectors 
k are chosen according to the periodicity condition 

2-JT 

k = —n , fieZ^ . (59) 
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The equal-time Green functions and the correlation dynamical equations of motion are 
expanded accordingly. 

In this section we investigate the behaviour of the system in the weak coupling 
limit by choosing certain initial conditions and propagating the connected equal-time 
Green functions according to the correlation dynamical equations of motion in different 
approximation schemes. 

The weak coupling limit is considered, since the properties of the regularized system 
with a given number of basis states are uniquely determined by the dimensionless 
quantity gL^^'^ (recall that in 2 + 1 dimensions has the dimension of energy); the 
limit gL^^"^ — > is therefore equivalent to the small volume limit. This should in turn 
make sure that the different numbers of basis states considered suffice to allow for an 
investigation of the asymptotic behaviour as the ultraviolet cutoff is taken to infinity. 

In order to investigate the nonabelian aspects we, furthermore, start with pertur- 
bative initial conditions that are close to the configuration given by the equal-time 
limit of ( P^^D at t = Iq. The advantage is that with (|T^), at t = to only the coupling 
terms contribute in the equations of motion, while the free (abelian) terms vanish (up 
to the explicit time dependence already present in the abelian ground state solution). 
Otherwise the time evolution in the weak coupling limit is usually dominated by the 



(uninteresting) propagation of a free field. The equal-time limit of (|l^) at t = to leads 
to 



2\k\ 



{At{k)A]{q))^ = 5'^%^.^ U.. - ^ 1 , (60) 



(n^(fc)n;(g-)). = 5^'5^,^.^ - ^ 



kj k'i \ \k\ 



(61) 



_,ab,, and 



(62) 



all higher order connected Green functions vanish. Note that while {AA)c, (nn)c, and 
{An.)c are purely transverse due to the gauge fixing condition (^) and the abelian 
limit of Gauss law (^), this is no longer the case for {IIA)^- 

While in principle the size L of the box introduces an infrared scale into the theory, 
( |60D - due to the periodic boundary conditions - still contains an infrared singularity 
which has to be regularized in some way. The task of finding gauge invariant infrared 
regulators for Yang-Mills theories is nontrivial and the resulting theory would in 
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general be hard to handle in the framework of correlation dynamics. We, therefore, 
resort to two different ways of handling this problem without having to introduce an 
additional mass scale by hand: we either leave out the zero momentum mode (which 
is not a gauge invariant prescription) or selfconsistently generate a mass for the initial 
condition by solving the quantum mechanical problem for the zero mode, only, which 
is described by the Hamiltonian 



H = in^n^ + -^e'""'e'"''A^A'jAfA'j , (64) 



with A1 = A1{k = 0), = Y\.1{k = 0) and [Af , U'^] = iS^^Sij. The Gaussian ansatz 

vi/(A)=n(^)'^'e-(^n^ (65) 

yields the configuration of lowest energy within the 2-point truncation scheme for 
a = 2 [l] hence an effective mass 

/ ^ \ 2/3 

mcff = 2" = ( ) ■ (66) 



Thus instead of (^) and (|6TD we use {uj{k) = yk^ + m% 



kikj \ 1 



{At{k)A]{q))^ = 5'^%^.^ 5,, - ) , (67) 



(nnfcln;(g-)), = 6^^%^., (5., _ ^ , (68) 

together with (^) and (^) as initial conditions and include = 0, such that at t = to 
the free propagation terms and the terms describing the coupling to the zero mode via 
the 4-point vertex cancel out. 

Before we begin our further investigation we note that, according to our experience, 
the solutions of the equations of n-point correlation dynamics with n > 3 (for any 
theory we have investigated so far, and also in the present case) become unstable at a 
certain point in time, when the absolute values of the equal-time functions suddenly 
start to rise steeply as a function of time. This does not seem to be a numerical 
effect, since the corresponding solutions remain unchanged if we use higher machine 
accuracies and smaller time step sizes. In the following we therefore focus on the time 
intervals where the solutions are still stable. 
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5.1 Infrared properties and gauge invar iance 

We now i) either leave out the zero momentum mode and initiahze (p(]|)-(p5|) or ii) 
include the zero momentum mode and initialize (0), (|68D, (|62|), and (^) at t = 
and then numerically integrate the correlation dynamical equations, which conserve 
the total energy in time. 

Fig. m shows the resulting electric field energy {ReEeiL, ImEeiL) 

Eel = {Eel) = \j d^x{Ii1{x)Iitm (69) 

as a function of time, where we have shifted the value at t = to zero. Since E^i 
is a gauge invariant operator, which is hermitean in the space of physical states, Ef.i 
should be a real quantity. However, the figures in the right column show that within 
all approximation schemes (first row: 2-point; second row: 3-point; third row: 4-point) 
we obtain an imaginary part that is roughly as large as the real part. These oscillations 
of Ef.1 in the complex plane are due to a violation of gauge invariance, leading to a 
non-cancellation of the terms in (p!7|) that differentiate between the Heisenberg and 
the Schrodinger picture in the propagation of gauge dependent quantities. While it is 
clear that within the 2-point truncation scheme (first row) gauge invariance must be 
violated due to the neglect of the 3-point vertex, our calculation shows that within 
the 3- and 4-point truncation schemes (second and third row), there is no significant 
improvement in this respect, which is partly due to a non- gauge invariant UV cutoff. 
One might argue that our initial conditions could also be a problem since they fulfill 
Gauss law only in the abelian limit, i.e. while at t = the Gauss law identities relating 
the 2- and 3-point functions are fulfilled, those relating the full 3- and 4-point functions 
are not. However, we have checked that by enforcing these identities at t = (cf. the 
discussion in context with Fig. we obtain almost no visible effect on our results. 

Some results concerning the infrared convergence are worth mentioning: Fig. |l| 
shows that within 2- and 4-point correlation dynamics (first and third row) the inclu- 
sion of the /c = mode (considering the 9 lowest lying momentum states) leads to a 
significant reduction of the amplitudes of the oscillation as compared to a neglect of 

— * — * 

the /c = mode (considering the remaining 8 lowest lying momentum states); in the 
4-point case the amplitudes are reduced by almost one order of magnitude. In contrast 
to that, within 3-point correlation dynamics the inclusion of the zero momentum mode 
only leads to a significant increase of fast oscillations with comparably small ampli- 
tudes, mainly in the real part of Ef.i, with frequencies given by the dynamics of the 
free field. The overall shape of the curves with respect to the slow, large amplitude 
oscillations - that violate gauge invariance - is almost unaffected by the inclusion or 
neglect of /c = 0. This can be viewed as an indication of infrared convergence within 
the 3-point truncation scheme, whereas there is no infrared convergence on the 2- and 
4-point level. 
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In Fig. 1^ we display the transverse and longitudinal contributions to the real part 
of Egi for the same calculation as in Fig. ^ the left column shows the results with 
the k = mode not included, the right column contains those with the k = mode. 
Within the 2-point truncation scheme (first row) we see that without the zero mode, 
the oscillations almost entirely (up to small deviations within linewidth) show up in 
the longitudinal sector of the theory. The inclusion of the zero mode then leads to 
oscillations in the transverse sector as well. This, however, is merely due to the fact 
that with our prescription for the longitudinal projector at k = 0, the zero mode is by 
definition completely transverse; the observed transverse oscillations are in fact entirely 
those of the zero mode. Within 3-point correlation dynamics the inclusion of the zero 
momentum mode leaves the longitudinal sector almost unaffected and only leads to 
an increase of fast oscillations already mentioned in the discussion of Fig. |l]. These 
oscillations apparently only take place in the transverse sector; a coupling between 
the transverse and the longitudinal sector only appears by generating an amplitude 
modulation for the fast, transverse oscillations; the transverse modes are of almost no 
effect on the longitudinal modes. Within 4-point correlation dynamics we observe a 
similar behaviour as in 2-point correlation dynamics, i.e. slow oscillations of the zero 
mode appearing in the transverse sector, with an additional increase of fast transverse 
oscillations as observed in the 3-point truncation scheme. 

The violation of gauge invariance in the present approximation schemes can be 
seen directly by looking at the violation of the Gauss law identities. The upper part of 
Fig. shows the quantity | Y^ki {^i{k)G"-{—k)) , which is a measure for the violation 
of the corresponding identity between 2-point functions of the form (AH) and 3-point 
functions of the form (AAII) (cf. eq. (0)). The calculations are the same as in Figs. 
^ and with the zero momentum mode included. Up to about t/L = 100, where 
the 3- and the 4-point approximation are still stable, all truncation schemes lead to a 
violation of roughly the same order of magnitude. In contrast to that, the quantity 

lEkM,,^,J^k,+k,+ks \e''''{A^{h)A'^{k2)G%h))\ displayed in the lower part of Fig. |, 
measuring the violation of the corresponding identity between 3-point functions of the 
form (AAn) and full 4-point functions of the form (AAAII) (cf. eq. (^)), becomes 
significantly larger with increasing order of the truncation scheme. Both results clearly 
indicate that the problem with gauge invariance cannot be cured by simply going to 
higher orders of correlation dynamics. 



5.2 Scaling behaviour and UV convergence 

In order to gain a better understanding of the behaviour of the theory within correlation 
dynamics, we now turn to an investigation of the system for different values of the 
dimensionless coupling gL^^"^ and for different numbers of plane wave states in the 
single particle basis. 

In the left column of Fig. ^ we again display the real part of the electric field energy 
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as a function of time, using the 21 lowest lying momentum states {k = included), 
for 3 different values of gL^^'^. In each of the 3 approximation schemes, the curves for 
the different values of the coupling apparently (up to the fast transverse oscillations 
in 3- and 4-point correlation dynamics, second and third row) have the same shape 
and are only squeezed in the direction of the time axis and stretched in the direction 
of the amplitude axis. This indicates that in the longitudinal (plus the zero mode) 
sector the theory exhibits a scaling behaviour with respect to the coupling. In fact, the 
corresponding frequencies (and amplitudes of dimension energy) behave like u oc g'^^^ 
as shown in the right column of Fig. ^ by rescaling all curves according to the scaling 
law 

/ \2/3 / \2/3 

E,i{g,,t) = i^) E,ii9iAf] t), (70) 

which is a property of any system with a purely quartic coupling. This can be seen as 
follows: consider a Hamiltonian of the form 

H = ln' + , (71) 

where this is meant as a purely formal notation with any indices or spacetime de- 
pendencies suppressed. One can now perform a canonical transformation according 
to 

A = g-'/^A' , n = g'/'W , [A, U] = [A', W] , (72) 
which leads to 

H = + \a'') (73) 

and thus to a g'^^^ scaling of all eigenvalues of the Hamiltonian. 

In case of 2-point correlation dynamics for Yang-Mills theory, a g"^^^ scaling can 
indeed be expected a priori for the longitudinal sector of the theory (if it decouples 
from the transverse sector) since the 3-point vertex is not taken into account and the 
2-point terms oc contain a transverse projection operator. The fact that the g'^^^ 
scaling survives in 3- and 4-point correlation dynamics despite the inclusion of the 



3-point vertex is, however, in line with the result obtained by Liischer that all 
eigenvalues of the (purely discrete) spectrum of Yang-Mills theory on a torus can be 
expanded in a power series of A = g'^^^. 

In order to investigate the asymptotic behaviour of the system as the UV cutoff is 
taken to infinity, the left column of Fig. |^ again displays the real part of the electric 
field energy as a function of time for different numbers of momentum states that are 
taken into account. In the first row we have plotted the corresponding result within the 
2-point truncation scheme and with the zero momentum mode excluded. All curves 
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apparently have the same shape and are only squeezed or stretched with respect to 
each other along the two axes. In the right column we show the curves to obey the 
relation 

Eei{n2,t) = ^^Y'E^iin,, t) (74) 

by rescaling them, where denotes the number of basis states taken into account; as 
in the case of Fig. |^ the rescaled curves are all on top of each other. This shows that 
the frequencies and the amplitudes of dimension energy diverge as the UV cutoff is 
removed - which is an unwanted behaviour since Yang- Mills theory in 2 + 1 spacetime 
dimensions is supposed to be finite. 

In the second row, the same results are plotted within the 2-point truncation scheme 
and with the k = mode now included. While the frequencies still scale like u oc n^^^ 
as in (ff^), the n^/^ scaling of the amplitude of E^i is destroyed by the inclusion of the 
zero mode. 

In the third row the corresponding results are shown in the 3-point approximation 
with the k = mode included; we have chosen only to display the longitudinal part 
in order to get rid of the fast oscillations in the transverse sector. From Figs. and |^ 
we already know that the results are essentially identical in this case if the zero mode 
is excluded. The rescaling of the curves works reasonably well in the first half of the 
displayed time interval; in the second half the application of ( fflD while still "going 
into the right direction" does not lead to a really impressive agreement of the rescaled 
curves. Furthermore, within the accuracy reached here, it only makes a small difference 
if in (^) the zero mode is counted or not. 

Relation (ff^) can be a priori expected within the 2-point approximation since there 
the longitudinal sector (if it decouples from the transverse sector) essentially behaves 
like theory with a purely quartic coupling, i.e. like a system described by the 
(momentum space) Hamiltonian 

H =lY.mm-k) + |^_E^ %+.^.+.^3+4,o'^(^i)<^>(^2)«f(^3)<f(^4) , (75) 

k kik2k-jk4 

where V denotes the quantization volume in any spacetime dimensionality. Computing 
the expectation value of ([75|) within a Hartree-Fock-Bogoliubov ansatz with effective 
single particle energies Q{k) then yields the energy functional 

W..EM.\Em)^ll^[E^^\ ,76) 



which is minimized by 

2V ^ n{k) 
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Eq. (|77D directly implies that Q{q) is independent of q and therefore leads to 

m^^ = 9'^'{E^ ' {^)"\ (78) 
which is the frequency scaling contained in ( [74|) and in ([70|). With 

E. = \^aCk) = \,^'4'L]\^f (79) 

\ k / 

we also recover the corresponding dependencies for the amplitudes. 

The result that (|7^ ) remains approximately valid in the 3-point truncation scheme 
indicates that, in contrast to the infrared behaviour, the cancellations required to guar- 
antee that the theory is UV convergent do not take place. This is again a consequence 
of the violation of gauge invariance. 

Furthermore, we note that within the 3-point approximation, in the equations of 
motion for the 3-point functions, the contribution of the terms coupling to a product 
of a 2- and a 3-point function via the 4-point vertex can to a good approximation be 
neglected. 



5.3 Restoration of Gauss law 



In Sections |5.1| and |5.2| we observed a violation of Gauss law, resulting in unphysical 
imaginary parts for gauge invariant hermitean quantities and in the non-cancellation 
of UV divergences. 

We therefore now try to improve our approximation scheme by partly replacing our 
dynamical equations of motion by the constraint equations arising from Gauss law (cf. 
Section There is of course no unique way to do this; we implement Gauss law by 
means of the replacement 

tk^u^ik) -^e'^'^ y: hA^,^Ai^imi<i2) , (80) 

91 92 

i.e. in an equal-time function that contains a longitudinal field momentum, the latter 
is commuted through to the right and the resulting n-point function is replaced by a 
sum over n + 1-point functions according to (|80|) (since we only use this prescription to 
replace 2- and 3-point functions, which possess no disconnected parts, we do not have 
to make the distiction between full and connected equal-time functions on the l.h.s. of 
(pOD). This is done throughout the time integration and the remaining functions are 
propagated as before using the correlation dynamical equations of motion. Of course 
( pOD cannot be applied for k = 0; the Gauss law for the zero mode, however, turns out 
to be fulfilled automatically within our investigation. 
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We apply this method within the 3-point truncation scheme, where we investigate 
the imphcations of a) either enforcing the Gauss law identities relating the 2- and the 
3-point functions (|43|) , (0) or b) enforcing the identities relating the 3- and the 4- 
point functions as well. In the latter case, the full 4-point function is given 

by products of two 2-point functions and we thus have a coupled system of equations 
for 2- and 3-point functions, which we solve by iteration. 

We mention that whenever there is a Green function containing more than one 
longitudinal field momentum, there is an ambiguity in the application of (^). In 
case of the 2-point/3-point identities, i.e. having a 2-point function containing two 
longitudinal momenta, this does not create any problem, since due to the symmetries 
of the system both ways for applying ( ]80| ) lead to the same result. In case of the 
3-point/4-point identities, this is not the case and we choose to symmetrize between 
the various possibilities; the stationary solution of our iteration scheme can then not 
exactly fulfill the corresponding identities. However, the remaining violation of the 3- 
point/4-point identities turns out to be very small compared to the original violation. 
Fig. P shows the electric field energy E^i and the magnetic field energy 
1 



E^ag = {E^ag) = ^J X {F^^{x) F^^{x)) (81) 

as a function of time obtained within the scheme above for 21 momentum states (A; = 
included). While the calculation that only implements the 2-point/3-point identity 
(upper part) becomes unstable at about t/L = 40, the main effect of the additional 
inclusion of the 3-point/4-point identity (lower part) is a stabilization of the system. 
In both calculations the real part of Emag on average decreases with time, while the 
real part of E^i on average stays constant (left column); both exhibit additional fast 
oscillations with frequencies that are essentially given by the dynamics of the free field. 
Since the total energy is given by Etot = Eei + Emag, this result implies that the price we 
have to pay for our restoration of Gauss law is the violation of energy conservation; the 
reason for this is the replacement of the 2-point function containing two longitudinal 
momenta, which contributes to E^i, by means of (|80|) . 

Energy conservation can in principle be restored by instead using the corresponding 
2-point/3-point identity to replace one of the 3-point functions (these do not show up 
in the total energy). Any corresponding prescription would, however, seem highly 
arbitrary, since there is a momentum summation and a contraction over vector indices 
on the r.h.s. of (^0|). Moreover, our attempts in that direction have lead to equations of 
motion that become unstable much earlier than the original ones and that still exhibit 
large imaginary parts for the electric and magnetic field energies. 

The imaginary parts of Eei and Emag (right column) both stay very small up to 
about t/L = 10, and then fast oscillations start to build up; if these are subtracted, the 
imaginary parts stay very small, which might be seen as an improvement. We note in 
passing that when neglecting the zero momentum mode, the fast free field oscillations 
almost completely vanish. 
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In Fig. 1^ the violation of the 3-point/4-point Gauss law identities is compared 
for three different approximation schemes with the same parameters as in Fig. |^: i) 
without implementing Gauss law, ii) implementing the 2-point/3-point identities and 
iii) implementing the 3-point/4-point identities as well. In the left picture we have 
plotted the quantity 1^^:,^^,^:^^., ^M.+k, e-''^{A-{k{)A]ik^)G\h)) 



Since there is 



no ambiguity with more than one longitudinal momentum, there is no such viola- 
tion if the 3-point/4-point identities are implemented. While the strength of the vi- 
olation oscillates without the implementation of Gauss law, it simply continues to 
rise if only the 2-point/3-point identity is implemented until the system becomes 
unstable at about t/L = 40. The same behaviour can be seen for the quantities 



6 ^ki,k2,kz,i,j ki+k2+k-i 



(middle picture) and 
(right picture), where the time inter- 



val has been chosen larger. One can see that without the implementation of Gauss law 
the equations of motion remain stable up to about t/L = 100. In the middle and the 
right picture, one can additionally observe the small violation arising in the scheme 
implementing the 3-point/4-point identities due to the above mentioned ambiguity in 
the apphcation of (^0]); the corresponding system of equations becomes unstable at 
about t/ L = 70. 

In Fig. 1^ we compare the time evolution of the real part of the magnetic field energy 
obtained with the implementation of 2-point/3-point and 3-point/4-point identities for 
different values of the coupling (upper part). After rescaling both axes (as in Fig. ^), 
all curves are on top of each other (lower part), indicating the validity of the g"^^^ 
scaling law ([701). The same result is obtained if only the 2-point/3-point identity is 
implemented. 

Fig. P shows the corresponding comparison for fixed coupling and different numbers 
of momentum states taken into account (upper part). After a rescaling of both axes (as 
in Fig. all curves are approximately on top of each other (lower part), indicating that 
the (I^fc)^''^ scaling law (74) holds as well. This result can also be reproduced within 



the scheme only implementing the 2-point/3-point identity. The ultraviolet behaviour 
is thus unchanged, the restoration of Gauss law has not helped to obtain the necessary 
cancellations and to achieve UV convergence. 

The behaviour observed in Figs. |^-|^ can again be reproduced by considering the 
simplified model for the 2-point dynamics of the longitudinal sector of the theory that 
is obtained from the ^''-like Hamiltonian (|7Sp. Within the 2-point truncation scheme 
the Gauss law identities (|4^) and ( ^3)) reduce to the abelian version of Gauss law, 
since the 3-point function is neglected; in our scheme given by (0) this corresponds to 
setting to zero all 2-point functions with longitudinal momenta commuted through to 
the right. In order to investigate the resulting behaviour within the $^-like model, we 
thus implement the conditions 

($(fc)n(-A^)) = iji{k)li{-k)) = , {U{k)<l>{-k)) = -i (82) 
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throughout the time evolution. The only remaining equation of motion then is 
dt{<l>{k)<^{-k)) = {U{k)<l>{-k)) + {<l>{k)U{-k)) = -i 

^ {^{k)^{-k)) = -It . (83) 



If we split up the expectation value of dTq ) into a part E^i corresponding to the elec- 
tric field energy and a part Emag corresponding to the magnetic field energy in the 
longitudinal sector of Yang-Mills theory according to 

E,, = i^(n(fc)n(-^)), 

k 

fclfe2fc3^4 

the insertion of (82) and (83) within the 2-point truncation scheme yields 

2 




E,i = 0, E_^ = -^ ^ . (85) 



Eq. ([85| ) is a particularly simple solution that fulfills the g"^^^ scaling law ( ffOD and 
the {Y^kY^^ scaling law (fflD , and that moreover yields no imaginary parts for E^i 
and Emag- Indeed, even though the 3-point functions are included, the magnetic field 
energies displayed in Fig. |]on average fulfill E^ag oc t^, while the electric field energies 
on average remain constant. 



In order to demonstrate the effect of the inclusion of the 3-point function. Fig. |TG 
shows the real part of the total energy obtained within the 2-point truncation scheme 
with the abelian Gauss law implemented in comparison with the total energy obtained 
within the 3-point truncation scheme with the 2-point /3-point and the 3-point /4-point 

— * — * 

Gauss law identities implemented. In case of the 2-point approximation, the k = 
mode is excluded, since otherwise the time evolution would be dominated by oscillations 
due to the zero mode dynamics. The inclusion of the 3-point function leads to a 
significant reduction of the violation of energy conservation; it, however, fails to lead 
to a complete cancellation of the corresponding terms. 



5.4 An extended momentum basis 

After having investigated the consequences of a restoration of the Gauss law identities, 
we now turn to a different method of improving 3-point correlation dynamics, which 
consists in altering the momentum basis. 

Within a finite momentum basis the violation of gauge invariance already shows 
up on the operator level, which means that the commutator of the Gauss law and the 
Hamiltonian becomes nonzero. This is due to the fact that for a number of terms 
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containing the projection operator in momentum space Y.k\k){k\ the completeness 
relation cannot be applied, since the sum of two or more momenta can lead out of the 
momentum space taken into account. 

By writing down the time derivative of one of the Gauss law identities relating 2- 
and 3-point functions, one can see that for these particular identities this problem can 
be cured by taking into account 3-point functions with one momentum outside the 
momentum space allowed for the 2-point functions. We, therefore, increase the allowed 
momentum space for the 3-point functions correspondingly. The remaining violation 
of the 2-point /3-point identities, which is then due to the 3-point truncation scheme, 
i.e. to the violation of higher order identities, is expected to be of a smaller effect on 
the time evolution than the original violation. 

We note in passing that it turns out not to have a significant effect if we also allow 3- 
point functions with 2 or 3 momenta larger than the 2-point cutoff, which is connected 
to the fact that the time derivatives of these 3-point functions then do not receive 
any contributions from the 3-point vertex term, which is the dominant nonabelian 
contribution (cf. the remark at the end of Section |5.1|) . 

We can also combine the theory with the increased momentum basis for the 3-point 
functions with the implementation of Gauss law identities discussed in Section p.3| . 

Fig. [TTl shows the electric field energy E^i (upper part) and the magnetic field energy 
Emag (lower part) as a function of time obtained within the above method by taking into 
account the 9 lowest lying momenta for the 2-point functions, which then (by summing 
up two of the momenta in a 3-point function) leads to an allowed space of 25 momenta 
for the third momentum (the sum of all three momenta has to vanish due to translation 
invariance). We compare the results obtained without an additional implementation of 
Gauss law identities and with the additional implementation of either the 2-point/3- 
point identity only or the 3-point/4-point identity as well, as in Section ^.3| . For both 
the real part (left column) and the imaginary part (right column) of E^ag we see that 
the results without Gauss law identities and those using the 2-point /3-point identity 
are almost identical. Up to about t/L = 10 the same holds for Eei. While for both 
approximations the system of equations becomes unstable at about t/L = 25, the 
additional implementation of the 3-point/4-point identities stabilizes the system up to 
about t/L = 80. 

Within all three approximations, on average, there is a large decrease in the real 
part of Eel and a comparably small decrease in the real part of Emag which implies 
that energy conservation is violated. The situation is similar to that of Section |5.3| , 
only with Eei and Emag exchanged. 

The imaginary parts of E^i and Emag stay very small until about t/L = 10. In 
case of no Gauss law constraints or only the 2-point /3-point identity implemented, 
the imaginary part of E^i then starts to build up fast (free field dynamics) oscillations 
with small amplitudes and in addition a slow large amplitude oscillation. If the 3- 
point/4-point identity is included, the slow oscillation vanishes. In the imaginary 
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part of Emagi the fast oscillations build up comparably large amplitudes, which are 
significantly reduced by the inclusion of the 3-point/4-point identity. 

In Fig. |l2]we demonstrate that the time evolution of the real part of the electric field 
energy again follows the g'^/'^ scaling of frequencies and amplitudes; the upper picture 
shows the corresponding curves within the approximation scheme implementing the 
2-point /3-point and the 3-point/4-point Gauss law identities for different couplings. 
The lower picture shows that the same curves with the axes rescaled according to ( [70D 
(as in Figs. § and ||) agree very well. The same result can be obtained within the other 



two approximations investigated in Fig. 11 



We note that at present we are not able to investigate the UV properties of this 
approximation since we cannot take into account enough basis states for the 3-point 
functions. 

Fig. compares the violation of Gauss law, measured by the same quantities as in 
Fig. ^ for all three approximation schemes discussed in this section and the unmodified 
3-point correlation dynamics discussed in Sections B?T] and ^71. In the beginning, the 
inclusion of higher momenta for the 3-point functions leads to a significantly smaller 
violation of the 2-point/3-point identity (upper picture) than in the unmodified case. 
The strength of the violation decreases, if more Gauss law identities are included, and 
there is a nonzero violation (becoming large as the system becomes unstable) even if 
the 2-point /3-point identity is implemented, which is due to the fact that the Gauss 
law for the zero mode, which is not enforced, is no longer fulfilled automatically as it 
is in case of Section |^ . 

The 3-point /4-point identity, if not enforced, exhibits a steeply rising violation 
(lower picture), which is much stronger than in the completely unmodified dynamics 
(Sections |5.1| and |5.2| ). If the 3-point/4-point identity is enforced, there is only a very 
small violation due to the zero mode. In view of these dramatic differences, it is 
remarkable how insensitive the results shown in Fig. |Tl| are to the inclusion of Gauss 
law identities especially in the early stages of the propagation. 



6 Adiabatic switching 

In Yang-Mills theory, as in any field theory, it is crucial to have a good approximation 
for the ground state of the theory. In this section, we therefore aim at obtaining an 
approximate solution for the equal-time Green functions in the ground state of SU(2) 
Yang-Mills theory on a torus in 2 + 1 spacetime dimensions within the framework of the 
equations of motion for connected equal-time Green functions as discussed in Section 

In our previous applications of n-point correlation dynamics to the nuclear many- 
body problem [^5|, ^ ^ and to $^-theory in + 1, 1 + 1 and 2 + 1 spacetime dimensions 
p8| , p9t pO| , we found that in the manifold of stationary solutions the total energy of the 
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system was in general not bounded from below, so that without additional boundary 
conditions the equations of motion cannot be used for a determination of ground state 
Green functions by a variational calculation. 

In the present case of equations of motion for Yang- Mills theory defined within the 
Heisenberg picture, the situation is even more complex, since in an equilibrium state 
only expectation values of gauge-invariant operators are stationary, and thus the task 
of finding the manifold of equilibrium solutions is significantly more involved than in 
one of the above non-gauge theories. 

In our previous investigations [JH, ?9, 28, |2^, 3^, in order to avoid the difficulties 



associated with a variational calculation, we successfully resorted to the method of 
adiabatically switching on either the coupling or the residual interaction, starting from 
the free solution or a stationary Hartree-Fock-Bogoliubov solution as initial condition, 
respectively. This method also has the advantage of requiring a significantly lower 
numerical effort than a variational calculation. 

However, in the following we will demonstrate that, at least within our approxi- 
mation scheme, the method of adiabatically generating a ground state configuration 
encounters severe problems in the present case of 2 + 1-dimensional SU(2) Yang-Mills 
theory on a torus. In the 2 + 1-dimensional continuum, without the presence of any 
regulators, it is obvious that the adiabatic method cannot work, since the coupling is 
the only dimensionful parameter and therefore the system must exhibit a discontinuous 
transition to the abelian limit. Even on a torus the theory can be shown to have such 
a discontinuous transition - in the framework of an implementation of a purely spatial 
axial gauge condition - due to the different numbers of independent zero modes for the 
Gauss law constraint in the abelian and the nonabelian case [^. Nevertheless, one 
could hope that in the presence of an additional infrared regulator (such as leaving out 
the zero momentum mode), this behaviour might change. Moreover, since the lowest 
energy state in the whole Hilbert space should automatically fulfill Gauss law [^, one 
could hope that in a possible adiabatic limit no additional difficulties arise by the fact 
that the propagation with a time dependent coupling, due to the induced explicit time 
dependence of the Gauss law operator, does not conserve gauge invariance. 

We proceed by parameterizing the coupling time dependently according to g{t)L^/'^ = 
at and start from the perturbative initial condition given by (|60|)-(|63|). We exclude the 
zero mode, since at t = we have ^'(t = 0) = and therefore the effective zero mode 
mass (^6]) vanishes, and we are left with the original infrared singularity of the free 
photon propagator. If the adiabatic method would work, then in the limit aL — * the 
system should propagate along an adiabatic trajectory. 

In order to investigate this we have plotted the real part of the total energy as a 
function of the coupling for different values of a in Fig. 0. In all three approximation 
schemes - unmodified 3-point correlation dynamics (first row), implementation of the 
2-point /3-point Gauss law identity (second row), and implementation of Gauss law up 
to the 3-point / 4-point identities (third row) - the curves, instead of converging towards 
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asymptotic values, apparently exhibit a scaling behaviour similar to the coupling or 
UV-cutoff scaling discussed in Section ^ (left column). The imderlying scaling relation 
can be shown to be 



Etot{(y2,g) = { —j Etot{oiu 



0:2/ 



3/5 



9) 



(86) 



by rescaling the curves accordingly (right column); the rescaled curves agree very well. 

Fig. ^ shows the same behaviour within the various approximation schemes taking 
into account an increased set of basis states for the 3-point function (cf. Section |5.4] ). 

Within (unmodified) 2- and 4-point correlation dynamics, the validity of ( ^61) can 
be seen as well, which is, however, not explicitly demonstrated here. 

All of the above results imply that the method of adiabatic switching does not work 
in the present case. 

In order to get an explanation for this behaviour we consider the simple example 
of a quantum mechanical system described by the (purely quartic) time dependent 
Hamiltonian 



X 



Ha{x,p,t) = y + («^)^y 



(87) 



Suppose we have found a solution of the time-dependent Schrodinger equation for 
« = 1: 

[^^t-H^{x,p,t)]^,{x,t) = 0, (88) 
and that we try to construct a corresponding solution for a 7^ 1 using the ansatz 
^^{x,t) = ^i{x',^t) (89) 
with x' defined by the canonical transformation 

x = /5-VV , p = (3^/^p'. (90) 
Inserting (p9| ) into the time dependent Schrodinger equations yields 
[idt - Haix,p,t)] ^„(x,t) 
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which directly implies a behaviour as in (^). If the initial condition is chosen inde- 
pendent of a, we furthermore have 



t = 0) = ^i{x,t = 0) = 0) ^„(a;, 0) = const (92) 

and thus our derivation only applies for the plane wave of momentum zero as initial 
condition, which is exactly what we have when starting with the unperturbed ground 
state. For plane waves with momenta ^ our consideration does not rule out (^), 
since the wavefunctions can always pick up phases that cancel out in expectation values. 

It is a priori clear that for the above example the method of adiabatic switching can 
never work, because a non-normalizable scattering state (e.g., a plane wave) cannot be 
continuously transformed into a normalizable bound state of a quartic potential. For 
this simple example, it is easy to circumvent this problem by simply introducing a mass 
(harmonic oscillator) term, which is switched off adiabatically as the quartic coupling 
is switched on. In the case of Yang-Mills theory, however, the situation is different: 
both states, the unperturbed one and the perturbed one, are non-normalizable due 
to their invariance under abelian or nonabelian gauge transformations, respectively. 
Therefore, the introduction of a mass term for the longitudinal Yang-Mills fields is no 
solution for the present problem since then the initial state would be a normalizable 
wave functional. 

7 Summary 

The main results of this work concerning the investigation of the weak coupling (small 
volume) limit of SU(2) Yang-Mills theory on a torus in 2 + 1 spacetime dimensions can 
be summarized as follows: 

• The violation of Gauss law within 2-, 3-, and 4-point correlation dynamics leads 
to large imaginary parts for the electric and magnetic field energies, which is 
clearly an unphysical result, whereas the total energy is conserved. 

• The frequencies of the slow oscillations in the longitudinal sector of the theory 
scale like g'^^^, as for a pure quartic oscillator. 

• Within 2- and 3-point correlation dynamics we also find the frequencies to scale 

/ \l/3 

like [J^kl , as in a $ -like theory described by a Hamiltonian of the form 

^IP + This result indicates that the cancellations necessary 

to ensure the ultraviolet finiteness of the theory do not occur on the 3-point level. 
This can also be attributed to the violation of gauge invariance. 

• Within the 3-point truncation scheme we find that the behaviour of the low 
frequency oscillations in the longitudinal sector is unaffected by the inclusion 
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or exclusion of the zero momentum mode, which is an indication of infrared 
convergence. 

• The modification of the 3-point truncation scheme by restoration of Gauss law 
identities of different orders does not produce any change in the coupling or 
ultraviolet scaling behaviour of the system; it, however, leads to a reduction of 
the unphysical imaginary parts of the electric and magnetic field energies, but 
also to the violation of energy conservation. 

• By taking into account an increased momentum basis for the 3-point functions 
we obtain similar results as for the restoration of Gauss law identities alone; the 
additional implementation of Gauss law identities, while stabilizing the system, 
does not lead to a significant change in the behaviour of the electric and magnetic 
field energies especially in the early stages of the propagation. 

• The adiabatic generation of a perturbed ground state from the free (abelian) 
ground state turns out to be impossible due to a scaling behaviour similar to 
that of a pure quartic oscillator. 

We thus conclude that standard many-body techniques such as the connected Green 
function approach in lower order (up to the 4-point level), which are successful for 
the nuclear many-body problem and for $^-theory, do not allow to compute the low- 
energy response of Yang-Mills systems, since either gauge invariance or the conservation 
of total energy is violated and since the UV divergences do not cancel out. This 
is mainly due to the fact that, for all approximation schemes studied in this work, 
the necessary cancellations between the 3-point vertex terms and the Hartree-Fock- 
Bogoliubov contribution to the 4-point- vertex terms (the dressed tadpole) do not occur. 



A A simple quantum mechanical example 

In this appendix we discuss a simple example, taken from ordinary quantum mechanics 
in 1 -|- 1 spacetime dimensions, in order to present the essential features of the approach 
of Section ^ to gauge fixing in the canonical formalism without having to deal with the 



complicated structure of Yang-Mills theory. The same example has been used in ^8 
to illustrate a different regularization procedure, which has the box normalization of 
plane wave states as quantum mechanical analogue. Both regularization schemes lead 
to the same result. 

We consider the Hamiltonian of a free nonrelativistic particle of unit mass (m = 1): 

"2 q2 

H=j, H^'^ix) = {x\H\^) = -y^(x) (93) 
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(in order to avoid confusion, throughout this appendix we denote operators by ■, their 
coordinate space representations are in analogy to Section ^ denoted by ■^). The 
eigenstates of (^Sj) are the non-normahzable plane wave states \p) with 



{x\p) = Ne'P'' . (94) 

We can now view p as the analogue of the Gauss law operator (^), i.e. we have 
\p,H] = 0. (95) 
The "gauge transformations" generated by p are the translations 
^(@) = e*ej5^ [/-(0)^(a;) = ^(x + e) , 

u~\e)xu{e) = x~e , u-\e)pu{e) = p . (96) 

If we define a space of "physical" states by the requirement of gauge (translation) 
invariance, the only state to survive is the plane wave with momentum zero, i.e. 

pl^p,,) = 0^\^pf,) = \p = 0) . (97) 

For the other plane wave states the momentum eigenvalue p can be regarded as the 
analogue to a static background charge distribution. 

Let Oinv be a translation invariant operator, i.e. [Oinv,p] = 0. In order to obtain the 
expectation value of Oinv in the physical state we first introduce the integral measure 
corresponding to (|]) by dividing out the volume of the translation group, 

djj,{x) = ( J dQ)^^dx , 

{^,H\d.n.\^,,) = J df^ix)%^ix)Ol,^,,ix) = ( J dQ)-' J dx Of„,l , (98) 

where the constant N is determined to be 1 by the requirement (1) = 1. However, 
dii{x) still leads to ill defined (infinite) expectation values of non-translation invariant 
operators, such as x^. In analogy to the derivation of ([13|) we thus write 



1 = A(x) j dQ6{f{x-Q)) , 
d^i{x)%y,{x)Otnv^ph{x) 

JdQj df,{x - e)A(x - Q)S [fix - 9)) %,ix - G)Or-®^,,(x - 9) 

dxA{x)5ifix))0tnvl (99) 



and define dv{x) = dxA{x)6 {f{x)), which in case of a function f{x) with only one 
simple root at X — Xq becomes 

0?Z/(x) = dx 6{x — Xq) . (100) 
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Instead of evaluating the spatial density corresponding to a translation invariant quan- 
tity by integrating over space and dividing by the volume, one thus can equivalently 
just pick the value at one specific point. 

The non-hermiticity of the Hamiltonian (^3]), when acting after a non-translation 
invariant operator within the regularization given by ( |100|) , can easily be discussed for 
the example of the operator x^: 

= J dx 6{x — xo) (^~2^x^ x^ = J dx 6{x — xq) (—1) = —1 

^ {H^p,\x'\%,) = J dv{x) {-]^dl%^{x)^ x^^!,h{x) = . (101) 

Furthermore, we can give a simple example to illustrate that only if x{t = 0) — = 
X — Xq appears to the left of all other operators (that it does not commute with), the 
gauge fixing condition = leads to vanishing expectation values: 

{^ph\{,x - XQ)px\^ph) = I dx 6{x - xo){x - Xo)- = 

f 1 1 

^ (,^ph\v[x - XQ)x\^ph) = J dx 5{x - xo)-{2x - Xq) = -Xq . (102) 

B Conserved colour structure of 4-point functions 

In this appendix we present the colour structure of the various connected four-point 
functions corresponding to a global colour singlet configuration, which is conserved in 
time by the equations of motion of 4-point correlation dynamics. We have 

{At^{xM1^{x,)Ar,{x,)Ai:{x4))c 

= d^^^'^d'^^'^^l - 5"i'^3)((A,,(fi)^..(^2)A.3(^3)^.4(^4)))r'* 
+5aia3^a,a4(i _ ^^i^^ ) ((A,^ (f , ) (f 3) A,, (xa) A,, (X4) 
^Sa,a,^a,as^-^ _ 5"^-^ ) (( A,, (f 1 )^M (^4) (^2) ^^3 (^^s) ))r'* 

+<5'^^'^^5'^^'^M'^^"M(Ai(xi)A,(x2)A,3(f3)A4(5^4)))r" , (103) 

mixMt^ixMt!{Ss)At:ix,)), 

= 6^^^'6''^'''{1 - 5'^^'^3)((n,,(xi)A,,(f2)A,3(f3)A,,(x4)))r'* 

+5'^^'^^5"^"M'^^"M(nil(^l)^.2(^2)A,3(f3)A4(^4)))r" , (104) 



/ c 

^aabb 
\aabb 
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+(5"^"^(5«^"M"^'^^((ni,(fi)ni,(f2)Ai3(:?3)A4(^4)))r" , (105) 



(n^;(fi)n^,^(f2)nr3=*(a;3)At'(^4))c 

= - 5«^-3)p^^(^^)n,,(f2)n,3(f3)Ai,(f4)))r'' 

- 5"^"^)((ni,(fi)n,3(f3)n,,(f2)A,,(f4)))r'* 

+5»^"^5«^"M"^"M(n,,(fi)n,,(f2)n,3(f3)^M(^4)))r" , (loe) 

and (n n n n)c has the same structure as {A A A A)c. As for the 2- and 3-point 
functions the double brackets {{■))c denote a function that does not depend on colour 
degrees of freedom. The superscripts ■"■"^"■^ .'^^''^ and ■"■^"■^ refer to the configuration of 
colour indices on the l.h.s. of the equation for the corresponding function to contribute, 
where the order of the colour indices refers to the order in which the corresponding 
operators appear in {{O O O 0))c on the r.h.s. (the notation becomes self-explaining 
by looking at the corresponding Kronecker-deltas in front of the colour independent 
functions) . 

In the case of a colour-singlet configuration, as investigated in the present article, 
one thus has 11 different 4-point functions that depend only on the spatial coordi- 
nates and on the vector indices of the corresponding field operators Ai and momenta 
Ilj. While the original, colour dependent connected equal-time 4-point functions are 
invariant under any permutation of operators (note that A^^x) and ^^j{y) commute 
in all connected Green functions of order > 3), the corresponding colour independent 
functions are invariant under a simultaneous exchange of two of the fields and the 
corresponding colour superscripts. 
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Figure Captions 

Fig.l : Electric field energy E^i = ^ Efc(nf (fc)n^(-A;)) as a function of time in 

2-point (row 1). 3-point (row 2), and 4-point (row 3) approximation with 
8 momentum states {k = mode not included) and 9 momentum states 
(A; = mode included); left column: ReE^i, right column: ImEgi 

Fig. 2: Real part of the transverse and longitudinal parts of the electric field en- 

ergy, El\ = i E,^ (5., - kAk)) (n«(^_)n«(-fc)), Ei^ = i j:ukAk)mk)U^(-k)), 

lij{k) = kikj/k!^{k: ^ 0), lij{k) = 0{k = 0), as a function of time in 2-point 
(row 1), 3-point (row 2), and 4-point (row 3) approximation; left column: 8 
momentum states {k = mode not included), right column: 9 momentum 
states {k = mode included) 



Fig.3: Violation of Gauss law identities measured by | • {A°;{k)G"'{—k)) 



[up- 



per part) and i Efc„fc„^3,,, %,+,^,+fc3 e^"^{At{k,)A'iik2)G%k^)) (lower part) 
as a function of time in 2-point, 3-point, and 4-point approximation with 9 
momentum states 

Fig.4: Real part of the electric field energy as a function of time in 2-point (row 
1), 3-point (row 2), and 4-point (row 3) approximation with 21 momentum 
states for different couplings (left column); curves rcscaled onto gL^^^ = 

0.01 via the assumption Eei{g2,t) = (^^^ Eei{gi, (^^^ t) (right column) 
Fig.5 : Real part of the electric field energy as a function of time in 2-point approx- 

— * — * 

imation without the k = mode (row 1), 2-point approximation including 
the k = mode (row 2), and in 3-point approximation including the k = 
mode (only transverse part of E^i, row 3) for different numbers n of mo- 
mentum states (left column); curves rescaled onto lowest number of states 

via the assumption Ef,i{n2,t) = (^^^ Eei{ni, t) (right column) 

Fig.6: Electric and magnetic field energy as a function of time in 3-point approxi- 
mation with 21 states using the Gauss law identities relating 2- and 3-point 
functions (upper part), and using the Gauss law identities relating 2- and 
3-point functions and 3- and 4-point functions (lower part); left column: 
real part, right column: imaginary part 



Fig. 7: Violation of Gauss law identities measured by 

6 '-^ki,k2,k3,i,j ki+k2+k3 



s-''{A-{k-,)A';{k2)G'ik,)) 
e'''''{U'}{h)A'{k2)G%h)) 



(left part), 
(middle part). 
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^^^\l:ua,S,.^JM.+k. e'^'^{Yi1{k,)Il]{h)G\h)) (right part) as a funo 



tion of time in 3-point approximation with 21 momentum states i) without 
using any Gauss law identities, ii) using the Gauss law identities relating 

2- and 3-point functions, and iii) using the Gauss law identities relating 2- 
and 3-point functions and 3- and 4-point functions 

Fig. 8: Real part of the magnetic field energy as a function of time in 3-point 
approximation with 21 states using the Gauss law identities relating 2- and 

3- point functions and 3- and 4-point functions for different couphngs (upper 
part); curves rescaled onto gL^^"^ — 0.01 via the assumption Emag{92,t) — 

(^f'E^,,{g„(ff't) (lower part) 

Fig.9: Real part of the magnetic field energy as a function of time in 3-point 
approximation using the Gauss law identities relating 2- and 3-point func- 
tions and 3- and 4-point functions for different numbers n of momen- 
tum states (upper part); curves rcscalcd onto 9 states via the assumption 

Emag{n2,t) = (^Y^^ Emag{ni, {^f^'^t) (lowcr part) 

Fig. 10: Real part of the total energy as a function of time in 2-point approxi- 
mation with 20 states {k = mode not included) using the abclian Gauss 
law identities and in 3-point approximation with 21 states (A; = mode 
included) using the Gauss law identities relating 2- and 3-point functions 
and 3- and 4-point functions 

Fig. 11: Electric field energy (upper part) and magnetic field energy (lower part) 
as a function of time in 3-point approximation with 9 momentum states for 
the 2-point functions and 25 momentum states for the 3-point functions i) 
without using any Gauss law identities, ii) using the Gauss law identities 
relating 2- and 3-point functions, and iii) using the Gauss law identities 
relating 2- and 3-point functions and 3- and 4-point functions; left column: 
real part, right column: imaginary part 

Fig. 12: Real part of the electric field energy as a function of time in 3-point 
approximation with 9 states for the 2-point functions and 25 states for 
the 3-point functions using the Gauss law identities relating 2- and 3- 
point functions and 3- and 4-point functions for different couplings (upper 
part); curves rescaled onto gL^^^ — 0.01 via the assumption Ejnag{g2,t) — 



91/ 



Emagigi, (g) t) (lowcr part) 



Fig.l3: Violation of Gauss law identities measured by 



(upper part) and i Eg,,g„,-3,,, %,+,-,+,-3 e'^'%A^{k,)A%k2)G%ks)) (lower part) 
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as a function of time in 3-point approximation with 9 states and in 3-point 
approximation witii 9 states for the 2-point functions and 25 states for the 
3-point functions i) without using any Gauss law identities, ii) using the 
Gauss law identities relating 2- and 3-point functions, and iii) using the 
Gauss law identities relating 2- and 3-point functions and 3- and 4-point 
functions 

Fig. 14: Real part of the total energy as a function of the time dependent coupling 
in 3-point approximation (row 1), in 3-point approximation using the Gauss 
law identities relating 2- and 3-point functions (row 2), and in 3-point ap- 
proximation using the Gauss law identities relating 2- and 3-point functions 
and 3- and 4-point functions (row 3) with 20 states (A; = mode not in- 
cluded) (left column); curves rescalcd onto aL — 0.01 via the assumption 

Etot(a2,9) = {^fEt,t{a,, {^fg) (right column) 

Fig. 15: Real part of the total energy as a function of the time dependent cou- 
pling in 3-point approximation (row 1), in 3-point approximation using the 
Gauss law identities relating 2- and 3-point functions (row 2), and in 3- 
point approximation using the Gauss law identities relating 2- and 3-point 
functions and 3- and 4-point functions (row 3) with 8 states for the 2- 
point functions and 24 states for the 3-point functions (A; = mode not 
included) (left column); curves rescaled onto aL — 0.01 via the assumption 

Etot{a2,g) = {^f'Etotiau {^fg) (right column) 
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